ITP-Budapest Report 648 



O 

(N 

> 
O 

in 



X 

5-H 



6-loop anomalous dimension of a single impurity operator 
from AdS/CFT and multiple zeta values 



Zoltan Bajnok a and Omar el Deeb 6 
November 16, 2010 



a Theoretical Physics Research Group of the Hungarian Academy of Sciences, 
b Institute for Theoretical Physics, Roland Eotvos University, 
H-1117 Pdzmdny s. 1/A, Budapest, Hungary 



Abstract 

O |' Anomalous dimension of the simplest nontrivial single impurity operator in the f3 = | de- 

formed theory is determined at six loops from the AdS/CFT correspondence. Liischer correction 
i-C ■ is evaluated at next-to-next-to-leading order (NNLO) in terms of multiple zeta values. The result 

can be simplified into the products of simple zeta functions and the same form of the correction is 
, expected for the Konishi operator at six loops, too. 

> 

. 1 Introduction 

AdS /CFT correspondence [1J relates the string energies on the AdSs x S 5 background to the anomalous 
dimensions of gauge invariant operators in maximally supersymmetric four dimensional SU(N) gauge 
theory. The correspondence is particularly useful in the large N limit when it can be described by a 
two dimensional integrable field theory, see [2\ for a general introduction and references therein. The 
advantage of the integrable two dimensional point of view lies in its non-perturbative nature, which 
nevertheless, can be used to calculate perturbative quantities, too. Indeed, the leading finite size 
correction of a two particle state can be described in terms of the asymptotical Bethe ansatz and the 
generalized Liischer formulas [3J, which, when expanded in the coupling constant, provides the exact 
perturbative anomalous dimension of the Konishi operator up to seven loops. The direct perturbative 
gauge theory calculations are very cumbersome and have been evaluated for the Konishi operator up to 
4 loops [U [5] only, where the wrapping part of the correction matches exactly the Liischer type finite 
size correction [3J . Thus the integrability based "string theory" techniques provide a way to go beyond 
the available perturbative calculations and collect information about higher order results. The simplest 
situation in this respect is the calculation of the anomalous dimension of a single impurity operator. 
Such operator is trivial in the M = 4 SYM but has non- vanishing anomalous dimension in orbifold or 
/3-deformed theories. The leading and next-to-leading wrapping correction have been calculated in [6J 
and in [7] and the aim of the present paper is to extend the calculation to next-to-next-to-leading order 
(NNLO). This order is important in the sense that it is not affected by double wrapping corrections 
thus the calculation based on the Liischer formula provides and exact answer. 

The /3-deformed theory is an M = 1 supersymmetric, exactly marginal deformation of the J\f = 4 
SYM and has been always the testing ground of the AdS/CFT duality. The dispersion relation of the 
excitations was calculated exactly in [8], which together with wrapping corrections [UHO] provides the 
exact anomalous dimension of a single impurity operator Tr(X Z L ). On the string theory side this 
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operator corresponds to a one particle state in the su{2) sector in finite volume. The leading finite 
size effect follows from the momentum quantization, or asymptotical Bethe Ansatz equation [11 J. As 
the volume decreases Liischer/wrapping type correction becomes important [51 [61 RL2"] and for an exact 
description they have to be summed up. An educated way to take all finite size correction into account 
would be to extend the Thermodynamic Bethe Ansatz of the Af = 4 theory [T3 l fT4 l \15 \ \W \ \T? \ [TB I fT9] 
to the /3-deformed case. The first step into this direction is done: in [20J the authors extended the 
functional Y system relations, while in |21) the scattering matrix of the theory together with its twisted 
boundary condition were identified. 

We are interested in the j3 -deformed theory at a particular value of /3, namely for the simplest 
nontrivial (3 — | and evaluate the Liischer correction at NNLO. This value of (3 is interesting as it 
correspond to a particular orbifold theory, too. As it was conjectured in |22| the anomalous dimension 
of the su(2) particle (Tr(XZj) in the ft = ^-deformed theory coincides with the anomalous dimension 
of the sl(2) particle (Tr(DZ Z)) of the orbifold theory. A clear description of the relation can be 
found in [23J where a TBA equation was also proposed to describe the single impurity excitation of 
the orbifold theory. The proposed infinite coupled integral equations in principle exactly describe 
the energy of the one particle state for any volume and coupling. However, as it was shown in [23J, 
following the work of [24, 25, 26], it reduces in the weak coupling limit to the Liischer correction and 
that is what we are going to evaluate now. We will focus on the structure of the result, especially how 
the various terms of the perturbation theory appear in terms of multiple zeta values. 

Multiple zeta values (MZVs) are the generalization of zeta functions: 



and are trivially related to the values of nested harmonic sums taken at infinity. They appear at 
various places of physics and mathematics including number theory, combinatorics, quantum field 
theory, statistical physics. For a recent review see [57] and references therein. What we would like to 
describe here is how the various MZVs appear at a given order of the calculation. We will see that 
at leading order we have simple zeta function, the next-to-leading order double MZVs will appear, 
while the novel result shows the appearance of triple MZVs at NNLO. Not only the depth of the MZVs 
increases with the orders but also their transcendentality degree. There are several nontrivial relations 
among the MZVs of a given transcendentality degree and one can always choose a good basis |28| . 
Any basis above transcendentality 8 contain MZVs which are irreducible in the sense that cannot be 
written as products of simple zeta functions. As the anomalous dimension of the simplest nontrivial 
operator at NLO had transcendentality 7 such irreducible MZVs did not appear [7j. At the NNLO the 
transcendentality is 9 so we expect the appearance of irreducible MZVs. Let us anticipate our result. 

We calculate the anomalous dimension of the operator Ti(XZ) in the & = \ deformed theory, or 
equivalently the anomalous dimension of the Tr(DZZ) operator in the Z 2 orbifold theory. On the 
string theory side they correspond to the energy of a one particle state in finite volume L = 1. It is 
a very special state as it has a vanishing rapidity u = 0, (p = it), which is not effected by finite size 
corrections. In contrast, the energy is shifted by vacuum polarization effects as 



The asymptotical Bethe Ansatz energy Eaba is simply the dispersion relation of a standing particle: 
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Jl>j2>->j»>0 




E = Eaba + AE 




1 + 8g 2 - 32g i + 256g 6 - 2560.g 8 + 28672 5 10 - 344064 5 12 + O (g 13 ) 



while AE corresponds to the wrapping interactions and has the expansion 



AE = AE 4 g 4 + AE 5 g 10 + AE 6 g 12 + ... 
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In [6] the LO correction was calculated 



AE 4 = 128(4C(3) - 5C(5)) 

The NLO correction turned out to be [7] : 

AE 5 = -128(12C(3) 2 + 32C(3) + 40C(5) - 105C(7)) 

In the main part of the paper we evaluate the NNLO correction and obtain 

AE 6 = -128(-192C(3) + H2C(3) 2 - 592C(5) - 440C(3)C(5) - 322((7) + 1701C(9)) 

which completes the anomalous dimension of the operator up to six loops. One can observe that the 
transcendentality is 9 and that no irreducible MZV appeared. This is due to an unexpected cancellation 
during the calculation. Let me note also, that the form of the correction matches exactly with other 
6 loop results like the ones in |29[ [23] . 

The rest of the paper is organized as follows: In the next Section we present the Liischer correction 
for a one particle state in the (3 — ^ deformed theory by recalling the available formulas from the 
literature. Section 3 contains the results of the calculation. As we used a novel method to sum up the 
various terms we present the four and five loop results in this fashion and show how we obtained the 
result at six loops. We pay particular attention how the various MZVs appear and describe how they 
can be simplified. Finally we conclude in section 4 and explain how our novel results could be used. 
The Appendix contains some technical information, especially on how we summed up nested harmonic 
sums. 



2 Calculation of the anomalous dimension from AdS/CFT 

As the anomalous dimension of the single impurity operator corresponds to a standing particle state, we 
have to calculate the standard Liischer correction. It describes how the one particle energy (dispersion 
relation) 

E = E(p = 7r) = Jl + I6g 2 sin 2 (|) = y/l + 16.g 2 = e 

is modified due to vacuum polarization effects. The p = tt momentum corresponds to vanishing rapidity 
u = which is protected by symmetry. Consequently the leading finite size correction is the energy 
correction which originates from virtual particles propagating around the circle and can be written as 

AE =-J2 I §-sTr(S%l(q,0))e-^ L + 0(g 16 ) (1) 

Here we sum up for all bound-states of charge Q of the mirror model, whose momenta are q, Sq\ 
describes how they scatter on the fundamental particle and e(q) denotes their mirror energy. We have 
to expand this expression in NNLO in g. The above form exactly describes the energy correction up to 
the order g 14 only since at the order g 16 double wrapping effects will contribute, too. Let us analyze 
the g dependence of the various terms. 



Exponential factor 

The mirror energy has the following parametrization 



Q ( q ) _ z (q,Q) . „±,„ n\-1 + i 



z + {q,Q) 
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Matrix part of the scattering matrix 

The contribution of the scattering matrix can be factored as 

BTtOSgfaO)) = S scalar (q, 0)sTr(S^ x ( 9 , 0)) 2 

where the super-trace of the matrix part contain the contributions of all polarizations. The various 
polarizations of the mirror bound-states can be labeled in the super-space formalism as (w^w^^ , 
w J 3 wf~ 2 ~ J 6162, w 3 3 wf~ 1 ~ : ' 63, lUgiu^ 1_J 64) see [3U1 |3J for the details. The super-trace in the /?- 
deformed theory evaluates in the previous basis as 

Q Q-2 Q-l Q-l 

sTr(citrix(^, *)) = E *) + E 5B2 ^ *) + * E - * E *) 

j=o i=o j=o j=o 

The S-matrix elements can be extracted from [3J- In calculating the corrections for the su(2) repre- 
sentative neither of the S-matrix contributions depends on j and they read explicitly as 

z — x+ 771 z — x+ (1 — x z ) x+ 771 \ 772 

SFl 3 {z,x)= Z —^-^ ■ SF2 3 {z,x)= Z+ - X l { y X+Z ~\ X -^ 
z~ — x + rji rj2 z — x+ (1 — x z ) x+ r]i 772 

The appearing string frame factors can be written as |^ = y^fr an d (ft) = fr- The x ± parameters 
depend on the momentum the usual way 



x± (p) = COt f ±l ( 1 + . 1 1 + 16.g 2 sin 2 



which in our case results in 



x ± = ±x = ±— 1 + e ; — = -1 
Ag x+ 



With these variables the super-trace of the matrix part takes a particularly simple form 



We have checked that the matrix part of the sl(2) representative of the orbifold model gives the same 
result. There, the deformation is such that the fermions do not contribute as the two undeformed 
fermionic S matrix elements are the same [SJ 13 1] . 



Scalar part of the scattering matrix 

The scalar part of the scattering matrix of a charge Q bound-state can be obtained by multiplying 
the scalar factors of its individual scattering constituents. The charge Q bound-state composed of 
elementary magnons as z — (21, . . . , zq), such that z~~ — and Zq = z + and the bound-state condition 
is also satisfied z 4 + = z^ +1 . Thus the full scalar factor as the product of the elementary scalar factors 
turns out to be [3J: 



S scalar (z,x) = T[S sl ^(z l ,x) =TT e -2M**,*) fL x (2) 
1 1 zj - x \ — : 
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In calculating the Liischer correction we have to evaluate this expression when z is in the mirror 
kinematics [\z~\ < 1, |z + | > 1). The analytical continuation has been carefully elaborated in [32J : 

S-LlaM*) = VQ^ x ) SSV{2) ^ x ) ( 3 ) 
A particularly good feature of the formula is that both expressions depend on only and not on the 
individual z^. Explicitly the su(2) scalar factors read as: 

qsu(2)y x _ ( z+ ~ X-)(z+ - X+) (1 - ^=Q(1 " ^t) 

while for the case (a^j > 1 following [32J we can write 



-i log £,3,1(2, a;) = $(z + ,.t + ) - $(z + ,a; ) - , + $(2 ,x ) 



i l0g 



We have to be careful as the conventions of |32] are different from ours. To turn into our conventions 
one has to replace z^ — > and —> x T . There are integral representations for $ and^P as 

. f dwi / dw 2 1 1 . r(l + ig(wi + w^ 1 - w 2 - w^ 1 )) 
$>{x u x 2 )=i<b — <i) — log- 



Cl 27ri / Cl 27ri ioi - xi w 2 - x 2 T(l - ig(wi + w x 1 - w 2 - u> 2 1 )) 

. . /" dwo 1 , r(l + iq(xi + xV 1 — w 2 — w^ 1 )) 
V(x u x 2 )=i 6 — £ log — — ^ K^- 

J Cl 27Ti w 2 - x 2 T(l - ig(xi + x 1 - w 2 - w 2 i )) 

where the integrations are for the unit circle. They are well-defined provided none of the Xi lies on the 
unit circle. 

As by now we have collected all the necessary formulas we turn to analyze their weak coupling 
expansions. 

Weak coupling expansion 

Our aim is to calculate the weak coupling expansion of AE for L = 1. In doing so we decompose the 
integrand of the Liischer correction (p} 



into a simpler rational part 

2 



PM)- 4x2(1 + 2 ^-- + )) 2 f*- 



((xz-) 2 ~ l)(x 2 - (z+) 2 ) \z+ 

which contains both the matrix part and the rational part of the scalar factor, and into the more 
complicated S part: 

i log £(t?,Q) = 2($(z+,a; + )-$(z + ,x")-$(z",a; + ) + $(z",x~)) 
-$(z + ,x + ) + $(z + ,x~) - $(z~,x + ) + ^(z" ,x~) 
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We expand these functions in g 2 as 

P(q, Q) = P s (q, Q)g 8 + P 10 (q, Q)g 10 + P 12 (q, Q)g 12 - 

£(<?, Q) = 1 + E a («, Q)g 2 + E 4 (g, Q).g 4 + . . . 
The expansion of the rational part is quite straightforward and we obtain 

4096Q 2 (-l + g 2 + Q 2 ) 2 



P*(q,Q) 



{q 2 + Q 2 )Hq 4 + (-1 + Q 2 ) 2 + V(i + Q 2 )) 

P 10 (q,Q) _ 8(7 9 4 + 3(Q 2 -l) 2 + 10g 2 (l + Q 2 )) 



Ps(q, Q) (q 2 + Q 2 )(q 4 + (-1 + Q 2 ) 2 + 2q 2 (l + Q 2 )) 

Pi 2 (q,Q) _ 2(27 + 268g 2 + 704g 4 ) 16(-3g 2 + 2 g 2 Q) 

P W (q,Q) ~ (3 + 28q 2 + 64q^(q 2 + Q 2 ) + (1 + 4q 2 )(l + q 2 - 2Q + Q 2 ) 

I6(3q 2 + 2q 2 Q) 8(-81g 2 - 117g 4 + I6q 6 + 27q 2 Q 2 + 16g 4 Q 2 ) 



(1 + V)(l + q 2 + 2Q + Q 2 ) (3 + 16g 2 )(3 + 10q 2 + 7q 4 - 6Q 2 + 10g 2 Q 2 + 3Q 4 ) 

In expanding the \E' and <& functions we use the same method we used in [7]. The expansion of the 
^(xi, £2) functions for \x 2 \ > 1 (string region) reads as follows 

*(si,x 2 ) = - — (^(l-ig(x 1 +x^ 1 )) + ^>(l+ig(x 1 +xj 1 ))) 

■ 2 

a 3 

+^(*2(1 - tff(a;i + a?! *)) + * 2 (1 + *s(a?i + x x + ... 

where ^ n (x) = (-p)™(log(r(x)) are the standard polygamma functions. If |xi| > 1 then $>(x x ,X2) 
starts at g 6 . In the opposite case using the identity $(£1,2:2) — $(0, X2) — ^(x^ 1 , x%), being valid if 
\x\\ 7^ 1, we can calculate the leading expansion of $ as 

<S>(0,x) = -( lE g-2C(3)g 3 + ...) 

x 

Using functional identities valid for integer Q we obtained 

v I n\ 16 Q afa ( Q~ i( l- 2 \ , a (Q + iq-2\\ 
= -72-^2 2 )+M 2 )) 



and 



E 4 feQ) = -E 2 (g,Q) + (g2 + Q2)2 -^TQ?( g '( 2 H 2 ( 2 W + 

32(, 1 (^^) + 5 1 (^±^)) + 8(5 3 (^f^) + . 3 (^±^)) + 2 C (3) 



where S n (x) are the analytical continuation of the harmonic sums S n (N) = X)fc=i P» • They are related 
to the polygamma functions atQ 

= (-l)" +1 n!(C(n + 1) - S n+1 (Q - 1) 
Let us use now these expressions to calculate the four, five and six loop wrapping corrections. 



^■For n = one has to replace £(1) with 7^. 
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3 Summary of the results 

Once we know the expansion of all functions we can systematically calculate the wrapping corrections 
by evaluating the expression: 

AE = - E I T~ [5 8p 8 + S W (Pw + i%£ 2 ) + g 12 (P 12 + P 10 S 2 + P 8 S 4 ) + . . .] 
Q=i J Z7T 

Let us analyze the energy correction order by order. 

4 loop contribution 

The first non- vanishing term is 

The integrand can be decomposed as 

P(a Q) = V R ^ {Q) | V R ^ iQ) I R *~ iQ) I R »~ {Q) I ^ 8 ~ +(Q) I R * +iQ) 
f^(q-iQV ^(<l + iQ) j q-iQ-i q + iQ-i q-iQ + i q + iQ + i 

This form is actually quite generic, that is we can always expand the rational part of the integrand as 



ei,e 2 ={±} j= 

As the function P 8 is real R^(Q)* = Rgj(Q) and #s~~(<3)* = #8~ + (<2)> #8~~(<3)* = Rs + (Q)- The 
various coefficients can be expressed in negative powers of Q and 2Q ± 1 as 

256 512 512 . 512a 512i 512* 

^84(<3) = T^T - ^ T\ + TrTTT - TT ! #83 (<2) = 



2 (20 -1) (20 + 1) ' 83W7 Q 3 (2Q-1) 2 T (2Q + 1) 2 

+ _ 640 512 512 512 512 512 

#82(0) - -754 + TyJ 



4 Q 2 (2Q-1) 3 (2Q-1) (2Q + 1) 3 (20 + 1) 
640i 512a 512a 512a 512a' 512a 



^ l(Q) 5 + Q 3 + (2Q - l) 4 (2Q -l) 2 (20 + i) 4 + (20 + 1) 2 
512a' 512a 512a' 512a 



^ 8 (Q) (20 + 1) 4 + (20 + 1) 2 ; Rs (0) (20 -I) 4 (20 + I) 2 
We can further decompose any function as 

5 p± 4 D± 

R± - R± X 7?± X R± • p± _ V" 8 J'Q fc . Z?± 8 J>( 2 Q ±:L ) fc 

^Sj ■- ^j.Q "T" rt 8j,2Q+l f rt 8j,2Q-l ' n 8j,Q ~ Qk ' /t 8j,(2Q±l) ~ (20 ± l) fe 

fe=l ^ fe=l V ^ ; 

(It can be extended an analogous way for any i?2fcj)- We calculate the integral by residues in which 
we take into account the poles on the upper half plane. They are located at q = iQ which we call 
the kinematical pole and at q = i(Q ± 1) which are the dynamical poles. The contributions of the 
dynamical poles are canceled when summed over the bound-state contributions (0)- Thus we have to 
keep only the pole at a'Q which result in the following sum 

00 00 
AE 4 = -Y / Res q=lQ P 8 (q,Q) = -J2 #8i(0) = 128(4C(3) - 5C(5)) 

Q=l Q=l 
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In summing up the Rsi,2Q±i terms here and later on we use the following trick: 

oo oo 

E^,2Q-l(3)+^ 2 Q +1 (Q)) = ^ 

Q=l Q=l 

This means that in the sum only the Q~ n type terms survive, which result in the sum X]q=i Q~ n an d 
contains single zeta functions reproducing the result of [6]. 



5 loop contribution 

The second order correction has the form 



A ^ = ~E / p-( P io{q,Q) + P S (q,Q)Ml,Q)) 

Similarly to the 4 loop result in the first term only the residue at iQ contributes and we obtain a single 
sum: P w has a similar decomposition as P s and only the J2q=i Q n = C( n ) terms survive which 
results in 

OO „ j oo 

-E / ^^ (g,g) = -E ii, ioi,Q(Q) = - 128 ( 48 c( 3 ) + 40 c(5)-io5C(7)) 

Q=l J Q=l 

In calculating the integral of the term Ps(q, Q)^ 2 (q, Q) we can use the fact that Pg is real and replace 
in £2 the expression Si(— 1 + V' q ) with its complex conjugate S\(— 1 + ®^ lq ): 



16Q , nCi (Q-iq-2 
q 2 ~TQ 



and take the real part of the integral. Interestingly the integral is real in this case and we have 

J ^(P 8 (<Z,Q)£ 2 (<Z,Q)) = / ^(P S (q,Q)± 2 (q,Q)) = AE 5Q + AE 5{2Q±1) 

The replacement has the advantage that S\( — 1 + does not have any pole on the upper half 

plane opposed to Si(— 1 + ^±H). Thus calculating the integral we can take the residues at iQ and 
i(Q ± 1) only. We use the same method in calculating the sums for the 2Q ± 1 type terms as we used 
for A4: we evaluate the Q — 1 term in the (2Q — 1) type terms, then we add the Q — > Q + 1 shifted 
terms to the 2Q + 1 type terms. These terms include the dynamical residues and the R^j(2Q±i) ^yP e 
kinematical residue and result in 

A£ 5 (2Q ± i) = "^V - 6) + 2048C(3) 

. From the remaining RsjQ type terms at the kinematical residue we got summands of the form 
SqLi Q n S m (Q — 1) = C( n J TO )j f° r n + m = 4 and n + m = 6. Explicitly summing them up we 
obtained 

AE 5Q = 512(8C(2) 2 - 4C(3) 2 - 11C(2)C(4) - 8C(2, 2) + C(2, 4) - 16C(3, 1) +4((3, 3) + 10C(4, 2) + 20C(5, 1)) 

We can express the multiple zeta values in terms of products of simple zeta functions either by the 
online high precision calculator of MZVs: EZFACE, or analytically, by using the algebraic relations 
between them [331 EH] and the data mine for MZVs: [5S]. Combining all contributions we arrive at the 
final 5 loop formula 

AE 5 = -128(12C(3) 2 + 32C(3) + 40C(5) - 105C(7)) 
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6 loop contribution 

The calculation of the six loop formula reduces to the evaluation of 

AE 6 = -jr f ^(Pi 2 (q,Q) + P w (q,Q)Mq,Q) + P$(<i,Q)Mq,Q)) 

The various terms have different structure in terms of multiple zeta values. The first term can be 
evaluated just by taking the residue at iQ and summing up the X)q=i Q n = C( n ) type single zeta 
function contributions of R^ 2 1 q(Q)' 

128(576C(3) + 480C(5) + 392C(7) - 1701C(9)) 

In evaluating the term PiqS 2 we use the trick to replace S\(— 1 + S+ig ) with its complex conjugate 
S\(— 1 + ®~% lq ) and calculate the integral of P10S2 by taking into account the contributions of the 
residues at iQ and i(Q±l) only. The results can be decomposed into the contributions of the dynamical 
poles together with the Rioj(2Q±i) type terms of the kinematical pole and into the most complicated 
contribution of the kinematical pole Riqjq terms. This latter one can be written in terms of MZVs as 
Sq=i Q n Sm{Q — 1) = C( n i m )i f° r H + m = 4, 6, 8. Explicitly summing up this term we obtain 

-1024(48C(2) 2 + (16 - 15C(4))C(4) + 8C(3)(C(3) - 6C(5)) + C(2)(32C(4) - 105C(6)) - 8£(6) 

-48C(2, 2) - 80C(3, 1) + 8C(2, 4) - 8C(3, 3) - 40C(4, 2) - 88C(5, 1) 
+15C(4, 4) + 3C(3, 5) + 45C(5, 3) + 105C(6, 2) + 210C(7, 1)) 

The basis of the multiple zeta values contain an irreducible depth two element which can be chosen 
to be £(3,5), see [28]. Interestingly, however, using the relations between the MZVs, [28], we could 
express the total contribution in terms of products of single zeta values as 

(-12007T 4 - 80tt 6 + 27tt 8 + 60750C(3)C(5)) 

675 

The other remaining terms give single zeta functions only. There are two terms 



, Q(1 + Q) , T , (n , cgg Q C Q(l + Q) , T , , n , 4096 2048 , 



and 

- 131072 frw*'(« + 5» < 5 » 

whose summation leads to alternating Euler sums. We do not write them here as these terms will 
cancel when we combine them with other terms coming from Ps^4- 

In calculating Ps^U, the trick of replacing S n (—1+ ^Ip 9 ) with its complex conjugate iS* n ( — 1 + ~% lq ) 
and taking the real part of the result can be applied. This works equally for S\{— 1 + Ss2 }g 1 (_ \ + 
®^ lq ), but not for S\(— 1 + ®'t l9 )Si(—l + ^~ tg ) which is manifestly real and needs a special care. For 
this reason, we separate £4 as 

S 4 = 5Re (E 4 ) + 64S!(-1 + ^±^)S 1 (-l + ^*?) 

where in S4 the replacement can be applied and it contains the rational part, too. The advantage of 
using E4 is that it does not have any pole on the upper half plane, so in calculating its contributions 
only the kinematical (iQ) an dynamical poles (i(Q ± 1)) have to be taken into account. In contrast, 
additionally to these poles, Si(— 1 + ^±M) has simple poles at q = i(2n + Q) with residue 2 for any 
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integers n. Taking into account the additional pole contributions we encounter a double sum of the 
form: 



^128P 8 (i(2n + Q),Q)Si(n + Q-l) 

Q=l n=l 

We first separate the contributions coming from R£. q , R^j 2 q±i an d Rf^i R% ■ I n calculating the 
R^j q part we use the result from the Appendix to arrive at: 

2048(6C(3)C(4) + 21C(4)C(2, 1) - C(3)C(2, 2) + 8((2, 3) - C(2, 5) - 4£(3)C(3, 1) + 16C(3, 2) 
-4C(3, 4) + 24((4, 1) + C(2)(8C(3) - 19C(5) - 24C(2, 1) + C(2, 3) + 4C(3, 2) + 12C(4, 1)) - 10C(4, 3) 
-20((5, 2) - 35C(6, 1) + 16C(2, 2, 1) - 2£(2, 4, 1) + 32C(3, 1, 1) - 8C(3, 3, 1) - 20C(4, 2, 1) - 40C(5, 1, 1) 

In calculating the contributions of 2 q±i we use the usual trick to shift the 2Q — 1 type terms and 
combine them with the 2Q + 1 type terms, but additionally we combine this step with an appropriate 
shift ri — >■ ri + 1, too. The resulting expression, apart form a rational part, contains ^(Q + 1) with a 
quite long prefactor which we do not write out explicitly here, as it will nicely combine together with 
the contributions of the kinematical poles. The contributions of iif^, Rg can be evaluated along the 
same lines and we obtain 

262144Q(1 + Q) 3 

(1 + 2Q) 4 ^ 2 +yj { ' 

Finally we have to calculate the contributions of the kinematical and dynamical poles. The dynamical 
poles are simple and result in poly gamma functions of the form ^,-(5 + Q) for i = 0,1. Combining 
their contributions with (|5l6p we obtain the simple result: 

4096(28C(3) - 31C(5)) 

The calculation of the contributions of the kinematical pole (q = iQ) is the most complicated. It 
slightly simplifies when we collect the remaining contributions from Q and from the term coming 
from the ^(Q + 1) type term of the residues at i(Q + 2n). The emerging sums, apart from rational 
expressions which can summed up easily with Mathematica, contain the products of harmonic numbers 
multiplied with rational functions. Decomposing them we have terms of the form 

00 

^^^(0-1)^(0-1) (7) 
0=1 w 

for m + 7i = 5, 6, 7 or 8 and alternating sums of the form 
00 / 1 1 \ 

for 77,777 and k < 5. The latter one can be treated by the usual trick of shifting Q in the first term, 
and evaluating the = 1 term separately. The difference of the resulting shifted harmonic sums can 
be simplified into a single harmonic sum, which then can be evaluated in terms of MZVs of depth two 
containing also individual irreducible elements like £(3, 5). In the expressions ([7]) we can use that 

S n (Q - l)S m (Q 1) = J2 i -L + J2 J- g JL + Sn+m (Q 1) 

j=l J fc=l j=l J k=l 

This will lead to multivariate zeta functions of the form £(fc, n, m), £(k, m, 77) and £(k, 777 + 77). Unfortu- 
nately the emerging expressions are not very illuminating and are rather lengthy. Magically, however, 
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when we combine all the contributions the irreducible £(3, 5) term disappears and we arrive at the nice 
result of the 6 loop correction which is: 

AE 6 = -128(-192C(3) + 112C(3) 2 - 592C(5) - 440C(3)C(5) - 322<C(7) + 1701C(9)) (8) 

In structure it is very similar to the result of other six loop calculations, like the NLO result for the 
twist three |29) or single impurity |23| operators. Let us finally mention, that we evaluated numerically 
the integral form of AEq and the result perfectly agreed with the exact analytical result ([5} . 

4 Conclusion 

In the paper we have evaluated Luscher's finite size energy formula for the simplest case in the 
AdS/CFT correspondence at NNLO and determined the anomalous dimension of a single impurity 
operator in the /3 = | deformed theory at six loops. The Luscher formula expresses the energy cor- 
rection in terms of the scattering matrix of the physical particle and the mirror bound-states. In 
calculating the weak coupling expansion of the correction we have to integrate the expanded scattering 
matrix elements for the mirror momentum and sum over the full bound-state spectrum QeZ. The 
expansion of the dressing phase in the mirror-physical kinematics at order n contains the (n — l) th 
power of the polygamma function. When integrated over the mirror momentum by residues derivatives 
of polygamma functions evaluated at integer values appear. These can be expressed in terms of nested 
harmonic sums which after summation over the bound-state spectrum result in multiple zeta values. 
At the n th order of the expansion MZVs of depth n appear, additionally, the degree of the MZV is 
increased by 2 from order to order. Irreducible MZVs, those which cannot be written in terms of 
products of zeta functions appear already at depth 2 and degree 8. Although during the calculation 
of the six loop anomalous dimension of the single impurity operator irreducible MZVs appeared, but 
in the final form they canceled out, leaving only products of zeta functions. 

We have no explanation why such a cancellation appeared and it would nice to evaluate the seven 
loop correction in this setting to see whether a similar effect takes place at the next order, too. It 
would be even more interesting to calculate the eight loop correction, since there a new type of double 
wrapping effect will appear. Another interesting problem is to evaluate these corrections using Laplace 
transformations [7], as in this case integrals of polylogarithms appear. Then using the theory of motives 
just as for scattering amplitudes in [35] might explain the form of the correction. 

The developed techniques can be extended to calculate the six loop anomalous dimension of the 
Konishi operator, which corresponds to a two particle state on the string side. We expect that the 
form of the wrapping correction will be the same, only the actual coefficients will change. This result 
would be very relevant to check to proposed excited TB A equations |17[ [18] . Indeed, the comparison 
of the infinite coupled TB A equation to the Luscher correction at five loops [24j [251 E3j confirmed only 
part of the TBA equations, while a six loop analysis would check the full set. 

The form of the correction we obtained might shed some light on how a possible finite system 
of NLIE type integral equations would look like as the weak coupling expansion of the anomalous 
dimension is very similar in structure to the expansion of the BES equation [36J. 
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A Summing up nested harmonic sums 

Here in this appendix we describe how we summed up nested harmonic sums. 
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At 6 loops we have triple and double sums and here we focus on the triple sums only. Consider a 
sum of the form 

oo 1 cc 1 
Q=l n=l 

We claim that the result is 

f(a, b) = C(a)C(M)+C(a)C(&+l) 

6-2 

+ ^(-l) J C(fe - j)C(o. i + 1) - (-!)" (C(a, & + 1) + C(a, &, 1)) 



where the Euler-Zagier sums are defined as 

C(ai,a 2 , . . . ,a n ) = ^2 



h>h>...>j n >o^ J - ■ ■ ■ J " 



The derivation goes as follows: 

First we use that Si(n + Q — 1) = Y^j=? 1 ^ an< ^ decompose the sum as 



oo oo 



Now exchanging the orders of the sums and using that 

£ IS ^ J - g (£ ?n-«' - *> - — (j + g i J J 

we can write 

oo 1 / oo Q-i \ 6-2 

E^E^ = £(-l) fc C(«, fc + W - fc ) - (- 1 ) 6 (C(«K(6 + 1) + C(a, 6, 1)) 
q=i ^ \n=l i=i / fc=o 

An easier manipulation gives 

oo oo 1 

E^Err?^ Si ( fl+ ^ 1 ) = c(<oc(M)-c(o+m)-c(m,i) 

Q=l V „=1 ^ + ^ 

In the 6 loop calculation terms with 6=1 appear. They are individually divergent but their appearing 
combination is finite. 
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